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It is now generally believed that galaxies were built up through 
^ I gravitational amplification of primordial fluctuations and the sub- 

sequent merging of smaller precursor structures. The stars of the 
^ structures that assembled to form the Milky Way should now make 

Q"^ ! up much or all of its bulge and halo, in which case one hopes to 

find "fossil" evidence for those precursor structures in the present 
Ph| distribution of halo stars. Confirmation that this process is con- 

O ' tinning came with the discovery of the Sagittarius dwarf galaxy, 

+i . which is being disrupted by the Milky Way, but direct evidence 

! that this process provided the bulk of the Milky Way's population 

of old stars has so far been lacking. Here we show that about ten 
Kj ; per cent of the metal-poor stars in the halo of the Milky Way, 

' outside the radius of the Sun's orbit, come from a single coherent 

■ structure that was disrupted during or soon after the Galaxy's for- 

mation. This object had a highly inclined orbit about the Milky 
Way at a maximum distance of ~ 16 kpc, and it probably resembled 
the Fornax and Sagittarius dwarf spheroidal galaxies. 

Early studies treated the formation of the Milky Way's spheroid as an 
isolated collapse, argued to have been either rapid and "monolithic"^, or in- 
homogeneous and slow compared to the motions of typical halo stars'^. A 
second dichotomy distinguished "dissipationless" galaxy formation, in which 
stars formed before collapse^, from "dissipative" models in which the col- 
lapsing material was mainly gaseous^. Aspects of these dichotomies remain 
as significant issues in current theories^, but they are typically rephrased as 
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questions about whether small units equilibrate and form stars before they 
arc incorporated into larger systems, and about whether they are completely 
disrupted after such incorporation. Stars from Galactic precursors should be 
visible today either as "satellite" galaxies, if disruption was inefficient, or as 
part of the stellar halo and bulge, if it was complete. 

Recent work examined the present-day distribution expected for the de- 
bris of a precursor which was disrupted during or soon after the formation of 
the Milky Way. Objects which could contribute substantially to the stellar 
halo near the Sun must have had relatively short orbital periods. Ten Gyr 
after disruption their stars should be spread evenly through a large volume, 
showing none of the trails characteristic of currently disrupting systems like 
Sagittarius^. In any relatively small region, such as the solar neighbourhood, 
their stars should be concentrated into a number of coherent "streams" in 
velocity space, each showing an internal velocity dispersion of only a few 
kms~^. Objects initially similar to the Fornax or Sagittarius dwarf galaxies 
should give rise to a few streams in the vicinity of the Sun. 

The high quality proper motions provided by the HIPPARCOS satellite 
allow us to construct accurate three-dimensional velocity distributions for 
almost complete samples of nearby halo stars. Drawing on two recent obser- 
vational studies^' we define a sample containing 97 metal deficient ([Fe/H] 
< — 1.6 dex) red giants and RR Lyrae within 1 kpc of the Sun and with the 
following properties: 

1. HIPPARCOS proper motions are available for 88 of them^^; for the re- 
maining stars there are ground-based measurements^^; in all cases accuracies 
of a few mas yr"-*^ are achieved. 

2. Radial velocities have been measured from the ground, with accuracies 
of the order of 10 kms~^. Metal abundances have been determined either 
spectroscopically or from suitable photometric calibrations^'^'^^'^^. 

3. Galibrations of absolute magnitude My against [Fe/H] for red giants^^'^^ 
and RR Lyrae^^, allow photometric parallaxes to be derived to an accuracy of 
roughly 20%. These are more accurate than the corresponding HIPPARCOS 
trigonometric parallaxes, but still remain the largest source of uncertainty in 
the derived tangential velocities. 

4. We estimate the completeness to be of the order of ~ 92%, based on 
the fact that there are eight known giants which satisfy our selection criteria 
but do not have measured proper motions. 

We look for substructure in our set of halo stars by studying the entropy 
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S of the sample, defined as: 
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where the sum is over the Ap elements of the partition, the i-th element 
contains Ni stars, and is the total number of stars. In the presence of 
substructure the measured entropy will be smaller than that of a smooth 
distribution, and will depend on the details of the partition; some partitions 
will enhance the signal, whereas others will smear it out. If there is no 
substructure then all partitions will yield similar S values, and no significant 
minimum value will be found. 

We implement this entropy test initially by partitioning velocity space 
into cubic cells 70 km s~^ on a side. This choice is a compromise. It leaves 
a large number of cells in the high velocity range empty, but in the regions 
containing most of the stars, there are at least a few stars per cell. 

It is necessary to quantify the significance of any observed low entropy 
value relative to the distribution expected in the absence of substructure. 
Here we do this by generating Monte Carlo realizations which test whether 
the kinematics of the sample are consistent with a multivariate gaussian 
distribution^^. We calculate entropies for 10000 Monte Carlo samples on the 
same partition as the real data; only for 5.6% do we find values of S smaller 
than observed. We have repeated this test for many partitions, finding a 
large number with probabilities as low or lower than this. In particular, for a 
partition with a 250 kms~^ bin in v^j,, and 25 kms~^ bins in vr and v^, {vr, v^j, 
and Vz are the velocity components in the radial, azimuthal and 2;-directions 
respectively), only 0.06% of Monte Carlo simulations have S smaller than 
observed. In general cubic cells yield lower significance levels, suggesting 
that the detected structure may be elongated along v^p. We conclude that 
a multivariate Gaussian does not properly describe the distribution of halo 
star velocities in the solar neighbourhood. 

At this point the main problem is to identify the structure which makes 
the observed data incompatible with a smooth velocity distribution. A com- 
parison of the three principal projections of the observed distribution to 
similar plots for our Monte Carlo samples reveals no obvious differences. To 
better identify streams we turn to the space of adiabatic invariants. Here 
clumping should be stronger, as all stars originating from the same progeni- 
tor have very similar integrals of motion, resulting in a superposition of the 
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corresponding streams. We focus on the plane defined by two components 
of the angular momentum: and = -J + J^, although is not fully 
conserved in an axisymmetric potential. In Fig. la we plot versus J_l for 
our sample. For comparison, Fig. lb gives a similar plot for one of our Monte 
Carlo samples. For Jj_ < 1000 kpc kms~^ and | J2I < 1000 kpc kms~^, the 
observed distribution appears fairly smooth. In this region we find stars 
with relatively low angular momentum and at all inclinations. In contrast, 
for J_L > 1000 kpc kms^ , there are a few stars moving on retrograde low in- 
clination orbits, an absence of stars on polar orbits, and an apparent "clump" 
on a prograde high inclination orbit. 

To determine the significance of this clumping, and to confirm it as the 
source of the signal detected by our entropy test, we compare the observed 
star counts in this plane to those for our Monte Carlo data sets. We count 
how many stars fall in each cell of a given partition of this angular momentum 
plane and compare it to the expected number in the Monte Carlo simulations. 
We say that the z-th cell has a significant overdensity if there is less than 1% 
probability of obtaining a count as large as the observed N'^ from a Poisson 

distribution with mean (iV*) = iVgT^ ^ AT], where N] is the count in the 

i-th cell in the j'-th simulation, and A^sim is the number of simulations. We 
repeat this test is for a series of regular partitions of p x g elements, with 
p, q ranging from 3 to 20, thus allowing a clear identification of the deviant 
regions. We find a very significant deviation in most partitions for cells with 
J_L ~ 2000 kpc kms~^ and 500 < Jj < 1500 kpc kms~^; the probabilities 
of the observed occupation numbers range from 0.03% to 0.98%, depending 
on the partition, and in some partitions more than one cell is significantly 
overdense. 

Given this apparently significant evidence for substructure in the local 
halo, we study what happens if we relax our metallicity and distance selection 
criteria. We proceed by including in our sample all red giants and RR Lyrae 
stars studied by Chiba and Yoshii^° with metallicities less than —1 dex and 
distances to the Sun of less than 2.5 kpc. This new sample contains 275 giant 
stars and adds 5 new stars to the most significant clump in our complete 
sample. Of the 13 members of the clump, 9 have [Fe/H] < —1.6, whereas 
the remaining 4 have ([Fe/H]) ~ —1.53 ± 0.12, indicating that they are also 
very metal-poor. These stars are distributed all over the sky with no obvious 
spatial structure. 
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In Fig. 2 we highlight the kinematic structure of the clump in the extended 
sample. The clump stars are distributed in two streams moving in opposite 
directions perpendicular to the Galactic Plane, with one possible outlier. 
This star has vr = 285 ± 21kms~^, and we exclude it because its energy 
is too large to be consistent with the energies of the other members of the 
clump. The velocity dispersions for the stream with negative (9 stars) are 
0",^ = 30 ± 17, (Tr = 105 ± 16, (T^ = 24 ± 28 in km s""*^, whereas for the stream 
with positive Vz (3 stars) these are cr,^ = 49 ± 22, = 13± 33, ci^ = 31 ± 28 
in km s^^. An elongation in the f^-direction is expected for streams close to 
their orbital pericentre (the closest distance to the Galactic Gentre; compare 
with other plots of simulated streams^). 

The orbit of the progenitor system is constrained by the observed posi- 
tions and velocities of the stars. The orbital radii at apocentre and pericentre 
are rapo ~ 16 kpc and rperi ~ 7 kpc, the maximum height above the plane is 
Zmax ~ 13 kpc, and the radial period is P ~ 0.4 Gyr, for a Galactic potential 
including a disk, bulge and dark halo^. We run numerical simulations of 
satellite disruption in this potential to estimate the initial properties of the 
progenitor. After 10 Gyr of evolution, we find that the observed properties 
of the streams detected can be matched by stellar systems similar to dwarf 
spheroidals with initial velocity dispersions a in the range 12 — 18 km s~^ and 
core radii R of 0.5 — 0.65 kpc. We also analysed whether the inclusion of 
an extended dark halo around the initial object would affect the structures 
observed and found very little effect. We derive the initial luminosity L from 
L = L7(/gi^"^ xCx f'"^), where L* = 350 L© is the total luminosity of the 
giants in the clump in our near-complete sample, fs^^^^ 0.13 is the ratio 
of the luminosity in giants with My and {B — V) in the range observed to 
the total luminosity of the system for an old metal-poor stellar population^^, 
C* = 0.92 is our estimated completeness, and /™ ~ 1.9 x 10~^ is the fraction 
of the initial satellite contained in a sphere of 1 kpc radius around the Sun as 
determined from our simulations. This gives L ~ 1.5 x 10^ L0, from which 
we can derive, using our previous estimates of the initial velocity dispersion 
and core radii, an average initial core mass-to- light ratio M/L ~ 3 — 10 Tq, 
where Tq is the mass-to-light ratio of the Sun. A progenitor system with 
these characteristics would be similar to Fornax. Moreover, the mean metal 
abundance of the stars is consistent with the derived luminosity, if the pro- 
genitor follows the known metallicity-luminosity relation of dwarf satellites 
in the Local Group^^. 

The precursor object was apparently on an eccentric orbit with relatively 
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large apocenter. Given that it contributes 7/97 of the local halo population, 
our simulations suggest that it should account for 12% of all metal-poor halo 
stars outside the solar circle. Figure 2 shows that there are few other halo 
stars on high angular momentum polar orbits in the solar neighbourhood, just 
the opposite of the observed kinematics of sateUites of the Milky Way^°. The 
absence of satellite galaxies on eccentric non-polar orbits argues that some 
dynamical process preferentially destroys such systems; their stars should 
then end up populating the stellar halo. As we have shown, the halo does 
indeed contain fossil streams with properties consistent with such disruption. 
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Figure 1: The distribution of nearby halo stars in the plane of angular mo- 
mentum components, vs. = \J + o^^^ ^^ar complete sample 
(a) and for one Monte Carlo realization (b). Our Monte Carlo data sets have 
the same number of stars and the same spatial distribution as the observed 
sample. The characteristic parameters of the multivariate Gaussian used to 
describe the kinematics are obtained by fitting to the observed mean values 
and variances after appropriate convolution with the observational errors. 
We then generate 10000 "observed" samples as follows. A velocity is drawn 
from the underlying multivariate Gaussian; it is transformed to a proper mo- 
tion and radial velocity (assuming the observed parallax and position on the 
sky) ; observational "errors" are added to the parallax, the radial velocity and 
the proper motion; these "observed" quantities are then transformed back 
to an "observed" velocity. Velocities are referred to the Galactic Centre; we 
adopt 8 kpc as the distance to the Galactic Centre and 220 kms~^ towards 
galactic longitude I — and galactic latitude fo = as the velocity of the 
Local Standard of Rest. 



Figure 2: The distribution of nearby halo stars in velocity space and in the 
Jz ^ Ji. plane. Data are shown for our original sample (filled circles) and 
for the extended sample of more metal-rich and more distant giants^'' (open 
circles). Candidates for our detected substructure are highlighted in grey: 
triangles indicate more metal-rich giant stars at distances > 1 kpc, diamonds 
more metal-rich giants at < 1 kpc, squares metal-poor giants at > 1 kpc, 
and circles metal-poor giants at < 1 kpc. 
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